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Abstract: The friction in plasmas consisting of two species with different temper- 
atures is discussed together with the consequent energy transfer. It is shown that the 
friction between the two species has no effect on the ion acoustic mode in a quasi-neutral 
plasma. Using the Poisson equation instead of the quasi-neutrality reveals the possibility 
for an instability driven by the collisional energy transfer. However, the different start- 
ing temperatures of the two species imply an evolving equilibrium. It is shown that the 
relaxation time of the equilibrium electron-ion plasma is, in fact, always shorter than the 
growth rate time, and the instability can thus never effectively take place. The results 
obtained here should contribute to the definite clarification of some contradictory results 
obtained in the past. 
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Plasmas both in the laboratory and in space are frequently in the state of partial 
thermodynamic equilibrium [l]-[3] (i.e. with an initial temperature disparity of the plasma 
constituents). Collisions in such plasmas will after some time eventually result in equal 
temperatures of the species, implying an evolving plasma. However, there exists a long 
standing controversy in the literature, which deals with the effects of this temperature 
disparity on the ion acoustic (lA) waves. 

In Ref. [1] it is claimed that the corresponding energy transfer may result in the 
instability of the acoustic mode at large wavelengths (within the quasi-neutrality limit), 
and that this growth may be described within the fluid theory. The necessary condition 
for the instability obtained in Ref. ^ for an electron- ion plasma is, in fact, very easily 
satisfied because it requires only a very small temperature difference between the two 
species (electrons and ions), viz. Te > 4Tj/3. This instability condition is obtained by 
using the energy equations including the source/sink terms originating from the collisional 
transfer, together with the corresponding friction force terms in the momentum equations. 
The sufficient instability condition is stronger because of additional dissipative effects, like 
viscosity and thermal conductivity. 

However, the current-less instability described in Ref. [4j is based on a model which 
disregards the same temperature disparity in the description of the equilibrium, which, 
due to the same reasons, must be time evolving. In other words, the effects of collisions 
in the equilibrium have been explicitly neglected. These effects have been discussed in 
Ref. [S], published one year after Ref. [1] and for the same quasi-neutrality case. There, 
it is claimed that there is no instability for any temperature ratio of the two plasma 
components, and moreover, that this holds even in a current-carrying plasma, as long as 
the difference between the electron and ion equilibrium velocity remains below the sound 
speed. All that was needed to come to that conclusion was to let the equilibrium plasma 
evolve freely in the presence of the given temperature difference. However, we observe 
that Ref. has apparently remained almost unnoticed by researchers, in contrast to the 
widely cited Ref. see e.g. Refs. [6]-[T0] and many others. 

In the present work, this controversy is revisited for any two-component plasma. Es- 
sential for the problem is the energy equation describing the temperature variation. In 
the simplified form that we shall use, it contains only the collisional energy transfer 
source/sink term on the right-hand side. This simplified form is used for clarity only 
because, according to Ref. 0], in the absence of currents, that term alone is supposed 
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to yield an instability. In view of the controversy mentioned above, here we give some 
details following Braginskii [11], where the energy equation for any species a is given in 
the form: 

2^a-Q^ + UaTaS/ ■ Va + -na{Va " V)Ta = Qa- (1) 

The corresponding equation for the species h has the same shape, but with a minus sign 
on the right-hand side. We use the Landau formula for the energy transfer source/sink 
term [12], Qa = 3mbUbanb{Tb - Ta)/ma, where [I3] 



= 4 ^''^ (Qaqb\ 2 riaLba . . 

''''' UJ Weo) 3{Tb + Tamb/mar/^' 

The Coulomb logarithm is given by Lba = \og[rd/bo\, = rdardb/{rla + ^Ibf^^ ^dj = 
v^j/upj, and Bq = [\qaqb\/{4:TTeo)]/[3(Ta + Tb)] is the impact parameter. 

The additional source/sink term of the form Ffa ■ [va — Vb), where Fja is the friction 
force acting on the species a, in the absence of equilibrium currents/drifts, is in fact 
nonlinear and will not be discussed here. The other (sink) terms, due to viscosity and 
thermal conductivity, are omitted only for the sake of clarity, i.e. in order to demonstrate 
more clearly the effect of the disputed collisional energy transfer term. The effect of 
these omitted terms is easily predictable. Equation ([1]) is valid for any species a, b, thus 
including the electron-ion plasma from Refs. fJ]-[9j. 

a) In Ref. ^ the collisions in the equilibrium were explicitly ignored. In that case, the 
two energy equations without the equilibrium effects, corresponding to the model from 
Ref. 0] read: 

^7 <~ o-'(a,fe)0V ■ ^(a,fe)l = ±^ ^ba U 61 — ^ al j 

(jZ O 

nib T^bi 

±2l'ba [TbG-Tao) . (3) 

nia no 
Here, the minus sign applies to the species h. 

The two momentum equations and the two continuity equations have standard forms 
and there is no need to write them down here. We stress only the presence of the friction 
force terms in the momentum equations. These are of the form Ffa = —maUaVabiva — Vb) 
and Ffb = -rribUbi^baivb - Va), respectively. 

In the case of quasi-neutral perturbations, the two number densities n(^a,b)i are calcu- 
lated from the continuity equations and are made equal assuming a quasi-neutral plasma, 
like in Refs. [1], [5] (this is typically done when dealing with wavelengths that are much 
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longer than the Debye length). The dispersion equation reads: 

iAmbiyba\ ( 2 5 a^lio + 



rua J \ 3 nia + mbJ 
Hence, even using the same model as in Ref. [4j, we conclude that there is neither an 
instability nor damping of the acoustic mode, regardless of the ratio Tao/Tt,o. 

Note that the momentum conservation condition = ^b^b^ba/ {f^ana) is nowhere 
used in the derivation of Eq. (jlj). This is because the friction terms vanish in any case. 
In fact, from the two continuity equations we have the velocities Vji = uriji/ {krio). These 
expressions, together with the assumption of quasi-neutrality, cancel the friction com- 
pletely. This remains so for any two species a and h as long as their charge numbers Za 
and Zb are constant. Further in the text we assume singly charged species. 

b) The derivations are now repeated for isothermal quasi-neutral perturbations. In 
addition, the energy equation may be omitted in the equilibrium also assuming that the 
relaxation time for the equilibrium temperature is much longer that the period of wave 
oscillations. Keeping the full friction force Fj in both momentum equations, and within 
the same quasi-neutrality limit, yields a rea/ dispersion equation oj'^ = k'^{Tao + Tbo)/{ma + 
nib). In the given limit the collisions (through friction) do not affect the isothermal ion 
acoustic mode. This fact is usually overlooked in the literature. The collisions appear in 
Eq. (jlj) only from the energy equations, yet they do not affect the lA mode. 

b.l) Using the Poisson equation instead of quasi-neutrality, for isothermal perturba- 
tions we obtain coupled and damped lA and Langmuir waves 



72/2 I 2\| 2 I 2" 
k \V^^ + j + ^pa + 

+ ^^^ra^Tfc + {^la^lb + ^Tb^pa) = 0- (5) 

In the collision-less limit the two modes ([5]) decouple by setting Ta = Tb, though strictly 
speaking this has not much sense because in this case the acoustic mode may lose its 
electrostatic nature, especially in pair-plasmas. For a pair (pair-ion, electron-positron) 
collision- less plasma the solutions are u"^ = cUp + k"^ (f^^ + f^j,) /2±[ujp + k^ {v^^ — v"^^) /4] . 
In the low frequency limit uj <^ ^p{a,b) and for an e-i plasma, from Eq. ([5]) we have 



^ — k'^v'l — i2Vf.iUjmf.r'^JP' / TTLi, so that the lA mode is damped 



u = ±kv, (l - rie"^] ' - iv,,Tl^k\ (6) 
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Figure 1: Approximative (full lines), and exact relaxation with time-dependent collision 
frequencies (dashed lines) of equihbrium temperatures ([7]). 



meZ/ei/mi, and 



Cg + f^j. The damping in Eq. 



is k- 



We have used 
dependent. 

c) From Eq. ([1]) it is seen that in a quasi-neutral homogeneous equilibrium, without 
flows/currents, the equilibrium temperature is also evolving in time as 

dT(^a,b)0 



r^, =1=2 — -J^baiTbo — Tao). 
at rria 



(7) 



Keeping the collision frequencies constant this gives the two temperatures T(^a,b)o = [^^(a,b)o(l+ 
exp(— 4z/a6t))+T'(a,6)o(l— exp(— 4z/abt))]/2 evolving towards the common value (Tao+Tl)o)/2. 
On the other hand, solving (JTj) numerically with time dependent collision frequencies ([2]) 
gives a slightly faster relaxation for the two temperatures. To get a feeling on the time 
scale, this is presented in Fig. 1 by taking uq = 10^^ m^^ and Tqo = 0.1 eV, T^q = STao- 

Eq. ([7]) is to be used in the linearization of Eq. ([T]) , which in the case Uao = = '^o 
yields: 



+ ?^«oV ■ Val = +2 — Pba {Tbl - Tax, 
at 6 VTia 

TTlb Hal - ribi 

rria no 

The corresponding equation for the component b is 

dTbi 2 ^ rrib , . 

^T- + • Vbl = -2 Uba [Tbl - Tal) 

at 6 m„ 



(8) 



(9) 



Here, in the process of linearization yielding Eq. ([9]), the term {3/2)nbidTbo/dt on the 
left-hand side, cancels out with the term —{mb/ma)i'ba(Tbo — Tao)nbi on the right-hand 
side after using the equilibrium equation ([7]) for the species b. 
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Hence, both Eqs. (jSj) and (Q are obtained taking into account the evolution of the 
equihbrium. There appears an additional asymmetry between the two energy equations 
(apart from the opposite signs of the first term on the right-hand side), due to the last 
term in Eq. ([8]). This extra asymmetry is a consequence of the fact that the internal 
energy of the two species may also change due to the presence of the new ingredient in 
the system, i.e., the perturbed electric field Uai — rihi = EqV ■ Ei/e (in the presence of 
the necessary collisions of course). However, it vanishes if the quasi-neutrality condition 
is used on the right-hand side in Eq. ([8]), which sometimes may be permissible in higher 
order terms but not in general, for example assuming that the source/sink term in the 
energy equations gives only small imaginary corrections to the frequency. 

However, regardless of the fact that the last term in Eq. ([H]) is used or not, the 
effects of the evolving equilibrium remain within Eq. ([H]) in both cases. Note also that 
the cancelation of the terms in the equation for the species b (which is due to evolving 
equilibrium as described above) remains intact. 

c.l) We stress that Eq. is obtained also by using Eqs. ([HI E]) in the quasi-neutral 
limit (implying that the last term in Eq. is omitted). Hence, the lA mode appears 
unaffected by friction in the quasi-neutral limit even if the energy equations are used and 
the equilibrium is described correctly as evolving. 

C.2) We now use the two energy equations ([HI IH]) with the Poisson equation. The 
dispersion equation becomes 
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Solving for the lA mode yields approximately the frequency of the lA mode 



(11) 
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The growth rate is: 
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2 ^2 + 2^ + (5/3)P« + vl,) - 2ul 



lA 



+2(1 - r„o/T,o) (z/a,A;2c^/a;,2^) [cu^^ - {5/3)evl, 



(12) 



In principle, Eq. (1121) reveals the possibility for a growing lA mode if the Poisson equation 
is used instead of the quasi-neutrality in a time-evolving plasma. For example, this can 
be easily demonstrated in the limit of negligible terms originating from the last term in 
Eq. dH]), i.e., on condition \(Tai — Tf,i)/{nai — nbi)\ ^ \Tao — Tf,o\/nQ, or in an alternative 
form, \(Tai — Tbi)/(Tao — Tbo)\ ^ r'^fjk'^\qb(pi\/Tbo. In that limit, the numerical solution of 
Eq. ffTOj) yields the growth-rate of the lA mode in an electron-ion plasma that is presented 
in Fig. 2. Here, no = 10^^ m~'^ and we take several values of T^/Ti, where Tj = 0.1 eV. The 
growth rate increases with T^/Ti but only up to T^/Ti ~ 3. For even higher values of the 
temperature ratio the instability ceases, this is represented by the dashed (Tg/Tj ^ 10) 
line. 

However, we stress that the system evolves in time, and in order to have a reasonable 
fast growth of the perturbations, the following condition must be satisfied [cf. Eq. ([7])] : 



Taking the electron-ion case like in Ref. |1] and the corresponding self-evident conditions 
nib <^ rua, Tao < Tbo, k'^v^^ < cu'^^ < uj^^ < cUpf^, from Eq. f[T^ to the leading order terms 
we obtain 



'Jr = 2{mb/ma)l^ba < 7. 



(13) 
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Figure 2: The growth rate of the lA mode in electron-proton plasma with no = 10^^ m 
and for several values of Te/Tj. 

Hence, because always ojpa > ^ia^ here we have 

7 < 7r-, (15) 

i.e., the system relaxes on a time scale that is (much) shorter than the eventual growth 
time, and consequently the assumed instability actually can not develop. We note that 
this is in agreement with some experiments, e.g. in a Q-machine plasma [Hj where the 
instability has never been observed even by cooling the ions to near room temperature 
while keeping various temperatures for the electrons. 

To summarize, the long existing controversy dealing with the stability of the ion acous- 
tic mode in plasmas in the state of partial thermodynamic equilibrium has been revisited. 
The results obtained here can be summarized as follows: i) The friction does not affect 
the lA mode in a quasi-neutral plasma; ii) Even using the non-evolving model equiva- 
lent to Ref. [1], there is no instability of the lA mode, contrary to claims from Ref. [1]; 

iii) When the equilibrium plasma is properly described as evolving in time, and as long as 
the quasi-neutrality is used, collisions do not produce a growht of the ion acoustic mode; 

iv) When the Poisson equation is used instead of quasi-neutrality, in principle there is 
a possibility for a positive growth-rate of the lA mode. It appears as a combined ef- 
fect of the breakdown of the charge neutrality from one side (introduced by the Poisson 
equation), and the heat transfer (the compressibility and advection in energy equation) 
from the other side, all within the background of a time-evolving plasma. However, as 
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the equilibrium plasma evolves in time, with the relaxation time given in Eq. ([7]), the 
obtained growth time must be (much) shorter than the relaxation time. Yet, this shows 
to be impossible and we conclude that there is no instability in the electron-ion plasma 
with an initial temperature disparity. 
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